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te LN RODUCT LON 


An investigation of the performance of a turbine stage 
when the rotor to stator axial spacing was varied ina sys- 
tematic way was carried out by Kane and reported in Reference 
1. The data on maximum efficiency vs. axial spacing obtained 
by Kane [Fig. 21, Ref. 1], was not sufficient to determine 
an optimum spacing because the minimum clearance was physically 
limited to 0.244 in. Therefore a recommendation was made in 
hoe. tO conduct further testing by modifying the rotor to 
allow a closer axial spacing. The original purpose in the 
current investigation was to conduct the recommended tests. 

A second purpose was to investigate ways to make the current 
rig/waterbrake system more stable. 

At the outset, significant revisions were required in the 
data acquisition programs to bring all data acquisition on- 
line for the first time. A first test was then conducted to 
checkout the rig and data acquisition system by attempting to 
duplicate run 7 of Ref. 1. Significant differences were 
noted in the interstage pressure ratios and the stator and 
rotor loss coefficients evaluated from the new data. In the 
attempt to resolve these differences, further tests were 
carried out which showed that the turbine stagnation tempera- 
ture had a measurable effect on the pressure levels (when 
properly non-dimensionalised) throughout the blading, whereas 


stagnation pressure had no measurable effect. It was also 
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evident in the early tests, that the speed fluctuations were 
much higher than the uncertainty interval (+1175 RPM) reported 
mo Ref. 1. 

The goal of the present study was to resolve the observed 
anomalies, and to obtain an understanding of the parameters 
affecting the accuracy of the data obtained from the rig. 
Three questions were addressed in the experimental program: 

(i) What was the status of the speed control using the 

current waterbrake and controller, and could it be 
improved? 

(ii) What was the cause of the observed effect of 

stagnation temperature? 

(iii) What was the effect of the uncertainties in the 
measurements on the accuracy and repeatability 
of data derived from the rig? 

The results of the investigation and responses to the 
above questions are contained in the main body of the report. 
Section II describes the test apparatus and associated instru- 
mentation, the data reduction procedures and test program, 
and gives an account of the procedures followed in each test 
run. Section III documents the results of the investigation 
which are discussed in Section IV. Conclusions and recommen- 
dations are contained in Sections V and VI, respectively. 

A detailed calculation of the uncertainties in parameters 
evaluated from the rig measurements is given in Appendix C. 
These calculations were important in that they illustrated 
clearly which measured parameters govern the accuracy of data 


obtained from the turbine test rig. 
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Ii. TEST APPARATUS AND PROCEDURES 


wae teot EQUIPMENT 

Figures 1, 2, and 3 show the test rig installation, Air to 
drive the turbine was provided by an Allis-Chalmers twelve 
stage axial compressor supplying pressure ratios of 3 to l. 
Pressure ratios of 4 to 1 were obtained (6 to 1 is available) 
across the turbine by sealing the turbine in the hood and 
actuating the exhauster assembly shown in Figure l. The 
exhauster assembly was an air ejector supplied from the Allis- 
Chalmers compressor. The exhauster reduced the pressure in 
the test hood to as low as half an atmosphere. Descriptions 
of the turbine test rig and associated instrumentation are 
given in Refs. 2-9. 

The geometry of the stage blading is contained in Ref. 
l. The different axial spacings possible with current 
hardware are shown in Figure 4. All tests in this investiga- 
tion were run in position 1. The waterbrake reported in 
Ref. 1 was modified twice during the course of the present 
investigation. First one "Slinger", and then a second, was 
added to attempt to attain more constant speeds. Figure 5a 
shows a side view of the dynamometer water wheel with two 
"slingers" installed. A close-up view of the water wheel 
with "Slinger" installed is shown in Figure 5b. Prior to 
installation of the slingers, the water level in the dynamometer 


Cavity would occasionally "collapse" inward, create a back 
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pressure, and cause a surge in RPM. The slingers were 


designed to prevent the water from "collapsing" inward. 


B. TEST INSTRUMENTATION 

A schematic of the system used to acquire the data is 
shown in Figure 6. The addition of the HG-78K scanivalve 
controller allowed the punched paper tape system used pre- 
viously, to be eliminated and provided complete on-line data 
acquisition and storage for the first time. A total of 60 
measurements were taken at each data point. Pressure measure- 
ments were taken using a 48 port Scanivalve. Fourteen other 
transducers were scanned for measurements of temperature, 
force, torque, and differential pressures. RPM was measured 
uSing an HP-5328A Universal counter. All measurements except 
RPM were made uSing an HP-3455A Digital Voltmeter. The counter 
and voltmeter readings were then channelled through an HP- 
3495A scanner to the HP-9830A computer and subsequently stored 
on a mass memory platter. The instruments and the computer 


were connected via the HP "interface bus" (HP-IB). 


C. DATA ACQUISITION AND REDUCTION 

Extensive modifications were required to the existing 
programs for data acquisition and reduction. A complete account 
of the modifications, descriptions of new software and cur- 
rent program listings are given in Appendix A. 

The following changes and additions were made: 

1) The data acquisition program was modified to incorporate 


On-line acquisition of Scanivalve and RPM data. 
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a.) 


ole) 


iv) 


v) 


The reduction programs developed by Solms [Ref. 7], 
Robbins [Ref. 8] and Boatright [Ref. 9] were modified. 
Software to sample, store and plot the RPM was 


generated. 


Software to determine the number of samples required 


to be taken of forces and torques was generated. 
Program "UNCERT" was written to compute the uncer- 
tainties in the interstage conditions and loss 


coefficients evaluated from measurements. 


D. TEST PROGRAM 


A list of the tests conducted is shown in Table I. Three 


series of tests were conducted, not always in sequential 


order, for the following purposes: 

1. To determine the number of samples required to 
produce a constant average in the data values from 
selected non-pressure channels. 

2. To evaluate RPM stability at various RPM's. 

3. To evaluate test rig performance at different RPM's, 
pressure ratios, and temperatures. 

The first test series was conducted at a P.R. = 3.0 and 11500 
RPM. The second and third series were conducted at a variety 


of RPM's and pressure ratios. 


At the beginning of each test, the pressure transducer 


in the Scanivalve and the differential transducer for the 


flow nozzle were calibrated using an air source and mercury 


monometer. The force and torque transducers were calibrated 
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by applying weights and adjusting the zero and range on the 
Signal conditioning units to give a digital output in 
engineering units. Using the HP9830A program "SCAMOD", 
given in Appendix A, data for a particular point on a given 
test was automatically scanned and stored. Average acquisi- 
tion time took 1 minute for a turbine performance point or 
65 sec. for RPM sampling. Each data point gathered on tur- 
bine performance was reduced immediately to determine if any 
obvious errors existed before taking the next point. At 
the end of the run, all turbine performance points were 
tabulated by batch processing. 

The performance tests were conducted holding stage pres- 
Sure ratio constant. This was usually done by holding PL 
at a desired level and varying hood (back) pressure to ee 


variations in the pressure ratio. However, in run 11, hood 


pressure was maintained at a constant level and PL was raised 
O 


to evaluate what effects, if any, changing stagnation pressure 
had on the results at constant pressure ratio. When tempera- 
ture effects were being evaluated, temperature was varied by 
closing off all, or a portion of, the water flow to the 
aftercooler and allowing the new temperature to stabilize. 


A change in total temperature, tT. , of no more than one degree 


O 
in five minutes was taken to indicate that conditions had 
Stabilizied. Tests to evaluate the RPM stability involved 
taking RPM samples at a variety of RPM settings. Each "sampling" 


consisted of 360 points gathered over a time interval of 65 


sec. The data were plotted and stored. 


20 





E. DATA STORAGE 

A list of the raw data parameters recorded for the per- 
formance tests and scaling factors applied to each is con- 
tained in Appendix B. Several channels were averaged after 
repetitive sampling to improve accuracy, and these are also 
noted in Appendix B. A list of tests conducted and files 
and record numbers where raw and reduced data are stored is 
given in Table II. Table A-XV contains a list of data reduc- 
tion plot routines by name, and the parameters which they 


mrot. 


F. DATA ANALYSIS 

Performance test data were routinely processed using the 
modified reduction programs given in Appendix A. 

Data from RPM sampling tests were plotted using software 
given in Appendix A. Also, too, RPM data sets, at different 
RPM's, were analyzed by Fourier analysis, using a Hewlett- 
Packard 9845A waveform analysis program. Each set of data 
was first transformed by removing the "DC" component, i.e., 
averaging 256 data points and then subtracting the average 
from each data point. The data were then fed, separately, 
into the program and, after transforming to the frequency 
domain, a power spectrum was plotted. The power spectrum 
used in the program employs Raleigh's Theorem of Conjugate 
Multiplication (see Ref. 10), to obtain the plot. After 
obtaining the power spectrum plots, the magnitudes of the 


first forty points were taken from each of the plots, 
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multiplied together, and used to produce a third plot. The 
third plot was designed to determine whether a natural fre- 
quency (and harmonics) were present at different speeds. 

A comparison of calculated and measured temperature 
drops was made using data from runs 10-15. The procedure 
consisted of taking the horsepower of the stage calculated 
from speed and torque measurements and computing the tempera- 
ture drop across the stage which would be required to produce 
that horsepower. The computed "delta T" was then compared 
to the measured "delta T" across the stage at the different 
stagnation pressures and temperatures of the tests. 

An investigation was conducted on runs 12 and 13 to 
determine what effects were evident on the reduced data of 
inserting a constant value of Wet = 1.02 lbm/sec. This was 
done by inserting a line 590 in program TTR2 setting Wor = 1.02, 
re-reducing the raw data from these runs and then batch 


processing the results. 
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IIL. RESULTS 


A. PRESENTATION OF RAW AND REDUCED DATA 

The results of the Fourier Analysis applied to RPM 
samples taken at 15000 RPM and 19000 RPM are given in Tables 
III through XII. The tables contain the raw data, data 
transformed by removing the "DC" component, time domain 
presentation, frequency domain presentation, and finally, 
the power spectrum data for each RPM. The correlation of the 
first 40 points of the power spectrum for the two RPM's is 
presented in Tables XIII and XIV. The raw correlation is 
contained in Table XIII, and the referred form in Table XIV. 
The referred form was obtained by dividing each value in 
Table XIII by the largest value (point 5). 

Table XV is a listing of the symbols used in Tables XVI 
through XXXVII. Figure 7 shows the pressure tap locations 
referred to in Tables XVII through XXX. The PaiLculeee and 
measured temperature drops from runs 10-15 are presented 
in Table XVI. The raw data for runs 10-15 are presented as 
voltages in Tables XVII through XXII. The performance data, 
engineering data, and converted raw data for selected parameters 
in runs 10-15 are presented in Tables XXIII through XXVIII. 
The results obtained by setting er 1.02 = constant are 
presented in Tables XXIX and XXX for runs 12 and 13. The 
results of tests to evaluate numbers of samples of key 


parameters to be averaged are presented in Table XXxXI. 
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B. PRESENTATION OF GRAPHS 

The graphs of RPM samples taken at speeds ranging from 

10600 to 20000 RPM are presented in Figures 8 through 20. 
The sampling rate is shown on each figure. Graphs of the 
power spectrum for 15000 and 19000 RPM are presented in 
Figures 21 and 22. The correlation of these two spectrums 
for the first 40 points of the frequency domain is presented 

in Figure 23. 

The graphs of four selected stage performance parameters 
uSing the data given in Tables XXIII through XXVIII are pre- 
sented in Figures 24 through 27. The effects of total tempera- 
ture and total pressure variations on Pi/SPL , stator nozzle 
taps, stator loss coefficient, and ee data parameters 
are presented in Figures 28-36. Analog records of the time 
variation in selected data channels are shown in Figs. 37-40. 

Finally, the various aspects of repeatability are graphically 
depicted in Figures 41 through 46 for runs 12 and 13 (both 


conducted at the same tT. . Py , and pressure ratio). 
O O 


C. UNCERTAINTY ANALYSIS 

An uncertainty analysis was conducted for run 10, point l 
to evaluate the uncertainty in stator and rotor loss coeffi- 
cients resulting from specified uncertainties in the measure- 
ments. The detailed analysis is presented in Appendix C. 
An outgrowth of the analysis was the discovery of the predomin- 
ant dependence of both loss coefficients on the accuracy of 


the calculated value of the resultant axial force. Asa 
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result, a simplification was found which allowed the calcu- 
lation of the uncertainty in the resultant axial force, 

Pi/P, pressure ratio, and stator and rotor loss coefficients 
to PON couneentca easily into a computer program, and made 

a part of the current reduction program sequence. Sample 
calculations showed only a slight deviation from exact results 
using this simplification (see Appendix C). The integration 
of this new program, called "UNCERT", into the reduction 
sequence is discussed in detail in Appendix A and Appendix C. 
It should be noted that the simplification used in the pro- 
gram underestimates the percentage error in stator loss 
coefficient by .03 times the actual error. In rotor loss 
coefficient the underestimate is .10 times the actual error. 
A complete exercise of this new program for runs 10-15 is 
presented in Tables XXXII to XXXVII. The uncertainties 
assigned to the raw data parameters presented in Appendix A 
were the same as those in Ref. 1. They were also consistent 
with the deviations foundin samples conducted during the 
present investigation. 


Of note in the results of the calculations was discovery 


Sem pegs ppm ct mms SSS Sea 
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coefficient. For the stator loss coefficient the uncertainty 
interval was approximately +.03 and for rotor loss coefficient, 
@eprtoximately +.05 to +.08. 

It was also discovered that the pressures on ports ll 
and 21 contributed collectively, fully 65% of the uncertainty 


interval in the resultant axial force. These channels 
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correspond to hood pressure and static pressure at stator 
tap "HUB #3" (or clearance-plate cavity pressure) respectively. 
An increase in accuracy and stability in these two measure- 
ments would improve significantly the uncertainty interval 
in both loss coefficients. It should be noted that the 
analysis in Appendix C takes account only of "normal" 
sampling uncertainty in the measurements and does not con- 
sider errors introduced through calibration, instabilities, 
Or repeatability problems. 

Finally, anomalies occurred in the raw data on runs 
13 and 15 which required altering the reduction sequence. 
The anomalies with modifications to the reduction program 


are listed in Table A-XxX. 
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IV. DISCUSSION 


A. RPM INSTABILITY 

The RPM stability was examined in early tests by 
digitizing repetitive samples using the HP Universal 
counter on-line to the HP 9830A and plotting the results. 
This capability, which had not existed previously, quickly 
showed that the RPM variations were unacceptable at all 
RPM's. Figure 8 illustrates the problem at a low sampling 
rate of 27 samples per minute. RPM deviations of greater 
than +/- 500 RPM were common. It was also found that 
changing the counter control setting in the data acquisi- 
tion program could drastically alter the appearance of the 
output. The longer was the sampling interval, the more 
stable appeared the RPM in the plotted output. Figures 8 
through 11 illustrate this effect. The data in the four 
figures were taken at the same RPM and only the counter 
setting was changed in the acquisition program. The graphs 
are labelled according to sample interval. There seems to 
be a relationship between the sample rate and the magnitude 
of the recorded variations. At a sample rate of 100 samples/ 
min., the variation was t 200 RPM. At a sample rate of 
6 samples/min., the variation was = .2 RPM. The most 
useful rate seemed to be 330 samples/min., which gave 
relatively good repeatibility (= 20 RPM), and gave a plot 
which qualitatively agreed well with independent observa- 


tions of RPM variations during the tests. 
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Figure 8 shows results with no "Slingers" installed 
on the waterbrake. Figures 9 through 12 were samples 
taken with one "Slinger" installed. The RPM control was 
still unacceptable in that RPM variations were far in excess 
of those required to obtain repeatable performance data. 
The addition of a second "Slinger" produced a marked change 
in stability. Figures 13 through 17 show the magnitude of 
the RPM excursions from the set point with two "slingers" 
installed. A comparison with Figures 10 and 12 illustrates 
the improvement obtained. The addition of the slingers 
clearly provided a substantial improvement in stability in 
the speed range up to 16500 - 17000 RPM. Between 17000 
and 18500 RPM, a range of severe instability still existed, 
as seen in Figure 17. The RPM instability seems to decrease 
at 19000 to 20000 RPM as shown by Figures 18 and 19, with 
20000 RPM being more stable, consistently, than 19000 RPM. 

These observations suggested that there might be a 
natural resonance encountered between 17000 and 19000 RPM. 
To evaluate this possibility a Fourier analysis was carried 
out on two sets of data taken at RPM's of 15000 and 19000, 
on each side of the instability. The power spectra cal- 
culated for the two sets of data are shown in Figures 21 
and 22. There were seen to be peaks in the low frequency 
end of the scale as was suspected from an examination of 
the RPM vs. time plots. The correlation of the two power 


Spectra was calculated anc is shown in Figure 23. A 
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definite peak at .107 Hz (period about 9.3 sec.) is evident 
with a second harmonic exactly at .214 Hz. It can be 
concluded that an unsteadiness at a frequency of .107 Hz 
exists which does not depend on the rotational speed. 
Furthermore, the low frequency of the instability suggests 
the problem lies within the controller or the waterbrake 
mechanics. Satisfactory speed stability is obtained up 

to and beyond the 17000 to 18500 RPM range. For the moment, 
it is recommended that this range be avoided when taking 
data. It is also recommended that further testing be 
undertaken to attempt to improve the speed stability in this 


range by small modifications of the waterbrake geometry. 


B. TOTAL TEMPERATURE AND TOTAL PRESSURE EFFECTS 

Runs 10 and 11 were conducted at different total 
pressures (59.7 and 74.5 in Hg, respectively), but at 
nearly equal values of the total temperature (635-639°R) 
and stage pressure ratio, (3.5). As can be seen in 
Figures 24 through 27, the various stage performance 
parameters, such as referred horsepower, efficiency, 
referred rotor torque, and referred flow rate, were 
measurably unaffected by total pressure changes. Also, 
referring to Figures 28 through 30, the stator nozzle 
taps and the stator exit pressure ratios are also unaffected 
by total pressure changes. It was concluded, therefore, 


that the total pressure had no measurable effect on the 
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interstage or stage performance parameters as long as 
stage pressure ratio remained constant. 

Runs 10 and Runs 12 through 15 were conducted at differ- 
ent total temperatures (574-676°R), but at nearly equal 
values of the total pressure (60 inches Hg) and the same 
pressure ratio (3.5). Figures 24 through 27 and 31 again 
show that there was no perceptible difference in overall 
stage performance. However, the interstage performance 
varies markedly. Figures 28 and 29 show that there was a 
difference in stator nozzle tap pressure ratios at taps 3 
and 4 in the supersonic portion of the nozzle. The locations 
of the taps are shown in Figure 7. Pressure tap 1 was not 
connected in these tests. Tap 2 is located at the throat 
of the nozzle, and taps 3 through 5 are located in the 
Supersonic portion of the nozzle. The stator exit pressure 
Cae PET’? 


Figure 30 shows that there was a definite change in each of 


Bacios, Puus/ Pt’ P/P , were also examined. 
the three pressure ratios with temperature. The results 
for two different stage pressure ratios, 3.5 and 4.0, 
respectively, are shown together with data from Ref. l 

in Figures 32 and 33. Runs 8 and 9 (TY =638°R) in the 
present test series differed from Runs 10 and 1l in Ref. l 


(Ty =670°R), in total temperature only. The same temperature 


O 

effect was evident at both pressure ratios. Finally, 
Figures 34 through 36 show the effect of temperature on the 
referred stator torque, resultant axial force, and stator 


loss coefficient, respectively. Again, the effect was 


clearly measurable. 
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In attempting to explain the variations in interstage 
parameters with total temperature, it should be noted that 
changes in stagnation pressure and temperature should 
produce no changes in any referred quantities (such as 
referred stator torque) at given values of pressure ratio 
and referred speed, unless geometry changes occur or the 
Reynolds number is particularly low. Also, at a constant 
total pressure, total temperature should have no effect on 
pressure ratios in the flow of a perfect gas through a 
De Laval nozzle. The observed pressure ratio differences 


in particular in the values of P,/P evident in Figures 


aa 
fees, and 32, in turn affect the Storrs tion of other 
interstage parameters such as the resultant axial force (RAF) 
and the. stator loss coefficient. In fact the differ- 

ences evident in these derived parameters with total 


temperature are directly related to the shifts in Pi/P, 
O 


with total temperature. Several avenues were investigated 

to explain the temperature effect. A summary follows: 
(l) Effect of Reynolds Number 

Differences in the nozzle boundary layer development 

resulting from the effect of stagnation temperature on 
the Reynolds number were considered. However, the 
Reynolds number variation obtained by changing stagnation 
pressure level more than overlapped the variation ob- 
tained by changing stagnation temperature. Since the 
results were unaffected by stagnation pressure changes, 


this could not be the mechanism. 
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(2) Non-Perfect Gas Effects 

The effect of temperature on the gas constant, 
ratio of specific heats and compressibility factor for 
air in the range of temperature of the tests, was found 
to have an insignificant effect on the predicted nozzle 
ratios. 
(3) Thermal Expansion Effects 

Geometry changes as a result of thermal effects 
could occur and affect the measurements in a variety of 
ways. It is important to note first that uniform expan- 
sion of the stator nozzles, so that A/A* remained con- 
stant, would leave the stator pressure ratio's unaffected. 

The observed decrease in static pressure ratio (or 
increase in Mach number) at increased stagnation tempera- 
ture, would imply that the throat area was required to 
expand proportionately less than than the downstream 
section of the nozzle. Calculations were carried out 
assuming no change in the throat size and linear growth 
of the aluminum nozzle material to the higher temperature. 
The change in area ratio was calculated to be less than 
0.003. At tap 4 in the stator nozzles, the pressure ratio 
macteased from 0.281 in Run 11 to 0.310 in Run 12. This 
corresponded to an area ratio decrease from 1.162 to 
1.120 (using isentropic gas tables), or an area ratio 
change of 0.042. This was an order of magnitude greater 


than could be explained by simple thermal expansion. 
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The fractional change in the throat area of the 
stator nozzles due to thermal expansion between the lowest 
and highest stagnation temperatures, would also be less 
than 0.003. An examination of the values of the referred 
flow rate, which is a measure of the throat area if the 
flow is choked, revealed that for all tests except Run 
12 they ranged within +0.005 of the value 1.025. fThis 
was consistent with the uncertainty calculated in Appen- 
dix C, resulting from the basic uncertainty in the indi- 
vidual measurements. In order to examine the effect of 
the small variations in ee on the calculated parameters, 
the data of Run 12 and Run 13 were re-reduced using 
. . = 1.02 in place of the measured flow rate (see 
Tables XXIX and XXX). A negligible effect on the values 
of interstage pressure ratio was observed. 

The effect of temperature on the labyrinth leak rate 
was questioned. In the existing reduction programs the 
effect of temperature on the "kinetic energy" factor 
had been included in the original calibration of Robbins 
(Ref. 8]. No correction was included however for the 
effect of temperature on the labyrinth clearance. When 
this was added (Appendix A), the effect on the results 
was found to be negligible. 

It was concluded that thermal expansion could not 


explain the observed effect of temperature on nozzle 


pressure ratios. 
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(4) Condensation 

It was questioned whether condensation could be 
occurring in flow in the supersonic portion of the nozzle. 
Condensation would certainly alter the static pressure 
levels, and condensation of the same air had been ob- 
served in the flow from a Mach 1.4 nozzle exhausting into 
the laboratory. A calculation was done to determine 
Static temperatures in the flow assuming an isentropic 
expansion of air and water vapor through the stator 
nozzle. The calculation is presented in detail in Appen- 
dix D. The calculation predicted supercooled conditions 
for the water vapor in the entire supersonic portion of 
the nozzle. Supercooling in excess of 100° F was pre- 


dicted at a stagnation temperature (Ty ) Oe al a Rema 


O 

: (670° R) resulted in supercooling of only 37° 
O 

F at the nozzle exit plane. Refs. 14 and 15 state that 


higher T 


supersonic flow containing water vapor, supercooled by 
90-110° F will result in condensation. The data for low 


Tr fall precisely in this range, so that condensation 


Eo 


might well have occurred. Inspection of columns P2 through 
P4 in Tables XXIII to XXVIII reveals that as total tempera- 
ture was raised, the values in these columns progressed 
non-linearly toward the values for the highest total 


temperature, (Run 15-T, = 676° R). 
O 


A test was conducted at a total temperature of 575° 
R and pressure ratio of 3.5 with the rotor removed, and 


hood exhausted to see if condensation could be observed 
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(visually) in the flow exiting the stator nozzle. No 

condensation was observed, but it should be noted that 
following the calculation in Appendix D, only 4.6% of 

the small amount of water in the flow must condense to 
effect observed changes in the stator nozzle pressure 

ratio at M= 1.2. The small scale of the nozzles, the 
swirling flow geometry and the background of metallic 

surfaces, did not help the detection process. 

A comparison of the temperature drops through the 
stage calculated from the torque-measurement of power, 
with the value measured using thermocouple probes is 
shown in Table XVI. It can be seen that the difference 
between the two measures is large at the low supply 
temperature, and small for the high temperature runs. 
Since the measurements of the referred power were in good 
agreement for all tests, the difference must be the result 
Sf an error in the outlet thermocouple indication, or 
in the assumption of perfect gas made in the calculation. 
Both effects would accompany the occurrence of condensa- 
tion. It was therefore concluded that the changes experi- 
enced in nozzle pressure ratios and stator exit pressure 
ratios were directly related to condensation effects in 
the nozzle. It is noted however that operation at the 
highest temperature gave consistent results, so that con- 
densation effects were absent, or minor. 

While testing at high total temperatures in Run 


15 (TY = 676° R), several pressure ports "blew" due to 
O 
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the plastic pressure tubing currently used; the plastic 
tubing softened, and failed at the elevated temperature. 
In order to obtain consistent, reliable results from the 
rig measurements, it is essential that high total tempera- 
tures be used to suppress condensation. It is therefore 
recommended that 1) the rig be reinstrumented using high 
temperature plastic tubing. (Plastic tubing is currently 
available commercially for use at temperatures up to 

250° F (710° R).) It is also recommended that 2) a 
minimum run total temperature no lower than 680° R, and 
preferably higher, be adopted for all performance testing 
in which interstage measurements are needed. To ensure 
that the highest total temperature is attained consistently, 
it is recommended that the aftercooler be drained when 


conducting such tests. 


C. CALCULATION OF THE REFERRED FLOW RATE AND LOSS COEFFICIENTS 
Sequential scanning of the fourteen nonpressure channels 

revealed that for a given RPM, temperature and pressure varia- 

tions were extremely small for normal RPM variations (+175 RPM). 

For this reason, the temperatures and all pressures except 

one were eliminated from the scan averaging techniques employed 

in Ref. 1 (see also Appendix A). It was found that channel 

24, the stator exit tip static pressure, was observed visually 

on the mercury manometer board to oscillate + 0.50 in Hg 

during tests. Channel 24 was therefore averaged when taking 


data (see Appendix A for full details). In addition, the 
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two force and two torque measurements were examined to 
determine how many samples per run point were required to 
produce a repeatable average value of the measurement. 

Table XXXI shows the results of this experiment. It can be 
seen that approximately 30 averaged samples produced a 
repeatable average in these parameters. On the basis of 
these results a 30 scan sequence was included in the acquisi- 
tion program. It is also of note that the variations in 
stator torque and resultant axial force after 30 scans are 
consistent with the uncertainty intervals assigned in Ref. l 
and in this report. 

Time traces of the analog output of axial force, closure 
force, stator torque, and rotor torque in comparison with an 
analog recording of the RPM are presented in Figures 37 to 
40. The traces were taken at 16000 RPM. The variations in 
both force measurements were seen to correlate well with RPM 
Variations and this implies that improvements in RPM stability 
will also result in reduced uncertainty in the computation 
of resultant axial force. It should be remembered however, 
that, as shown in Appendix D, the two pressure port measure- 
ments (ll and 21) overwhelmingly dominate in producing the 
uncertainty in the resultant axial force. The variation of 
stator torque with RPM was small and this was consistent 
with the 100 sample survey in Table XXXI. Finally, the varia- 
tions in the rotor torque were larger than in the stator torque 
which again was consistent with the results of Table XXxXI. 
However, the rotor torque variations were seen to correlate 


poorly with RPM variations and exhibited the most noise and scatter. 
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Three key parameters, and Cpr are calculated 


Wref! st 
in the course of the data reduction, and the accuracy obtained 
depends on many measurements. What follows is a discussion 
of the various factors affecting the calculation of these 
three parameters: 

A comparison of Runs 12 and 13 revealed that there were 


differences in both W and Cc which could not be readily 


ef 


explained by the differences in total temperature, total 
pressure, or by the possible occurrence of condensation. Both 
runs were conducted at identical temperatures (574° R), 
pressures (60 in Hg), pressure ratios (3.5), and relative 
humidities (70%). Inspection of Tables XXV and XXVI revealed 


approximately a 1/2 % change in W In addition, Figure 41 


reft° 


showed clearly that there were differences in See Calcula- 


=eloms of Je involve three measurements; the pressure and 


ef 


the temperature upstream of and the pressure difference 
@emess the flow nozzle, which is located in the line from 
the Allis-Chalmers compressor to the rig inlet plenum. By 
interrogating the reduction programs it was revealed that the 
effects of small (realistic) changes in nozzle pressure or 


bemperature affected the values of We very slightly. How- 


ef 


ever, the value of W was affected significantly by small 


ef 


(but realistic) changes in the pressure drop across the flow 
mee2gle. Point 3 on Run 12 and point 5 on Run 13 were com- 
pared. It was found that the AP across the nozzle was 20.884 
an HO for @he first pointeand 21.079 for the second point. 


This is a change of .195 in H.,O representing a .9% change in 


2 
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meme; but thas difference is sufficient to explain most of 


the difference in the W Inspection of all values of We 


ef 
femerkuns LO-15 revealed only + .002 lbm/sec variation in the 


ref~* 


magnitude in any given run. However, the average We changed 


ef 
slightly from run to run. Based on the comparison of Runs 

12 and 13, this would appear to be due to slight drifts in 
transducer output during different tests, or to minor calibra- 
tion errors either before or during the runs. It is there- 
fore recommended that the calibration of the AP nozzle trans- 
ducer be checked rigorously by the investigator at intervals 
Emeoughnout each test. 

A re-reduction of data from Runs 12 and 13 was carried 
out with ee = 1.02 lbm/sec = constant entered into the pro- 
gram as discussed in Section IV.B. This re-reduced data is 
presented in Tables XXIX and XXX. The original stator loss 
coefficient data are shown in Figure 4l and the re-reduced 
data are shown in Figure 42. It can be seen that, while the 
calculated stator loss coefficients were in better agreement 
than before, they were still higher in Run 13. The effects 
of Wot then do not account totally for the differences seen 
in oc (and CR): The data in Figure 43 illustrates that, 
although there were differences in Sc and - that again, 
stage performance parameters were nearly the same. Figures 


44, 45 and 46, containing data from Run 12 and Run 13 show 


that there was only a slight difference in measured stator 


Saye ee 


Seems ecm rem ce 
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and the calculated interstage pressure ratio (P,/Pt je 

The distribution of the data in Figures 45 and 46 se entirely 
similar. This result is in complete agreement with the 
results of the uncertainty analysis given in Appendix D 

which emphasizes the dominance of RAF on the calculation of 


P,/P, . This also points to the necessity for an extremely 


eeu calibration of the force capsules before each run. 

The stator hub and tip pressures, on ports 21 and 24, 
were observed (on the mercury manometer board) to exhibit 
unsteady behavior. The unsteadiness on Port 24 was qualita- 
tively worse than on port 21 by about a factor of two. The 
acquisition program was altered initially to average both 
these ports over 30°scans. However, this resulted in the 
Scanivalve having to cycle completely through 48 channels 
for each of the 30 samples. This procedure required 300 
cycles of the Scanivalve to obtain ten data points in each 
performance test. Excessive scanivalve wear would result 
from such a procedure. It was therefore decided to only 
average port 24 since it was observed to be the most unsteady. 
However, the uncertainty analysis revealed that port 21 exer- 
cises a far greater influence on the uncertainty in the resul- 
tant axial force than does port 24, by a factor of nearly 26. 
In selecting between ae two measurements, port 21 should 
have been averaged, not port 24. It is therefore recommended 
that a pneumatic damper be inserted in the pressure lines to 
each of ports 21 and 24. An improvement in the steadiness of 


the pressure at these two ports, particularly 21, should result 
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in a significant improvement in the steadiness and reliability 
of the calculated values of resultant axial force, interstage 
pressure ratio, stator and rotor loss coefficients. The 
proper averaging of the pressures on these two ports would 

do more to improve the reliability and smoothness of the 
output calculations than any other recommendation which has 
been made. Hood pressure, the other key input affecting 

the resultant axial force calculation, was observed to be 

very steady throughout all tests. 

In summary, it is recommended that, prior to each run, 
the investigator, with the cooperation of the test engineer, 
carry out a systematic verification of the calibration of 
the measurements shown to have a controlling effect on the 
reduced data. Where practical, calibrations should be 
verified at intervals throughout the test period. The key 
parameters and recommended procedures to follow are given 


in Appendix E. 
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V. CONCLUSIONS 


1. The speed control of the waterbrake and turbine 
system manifests a characteristic frequency of .107 Hz 
at peeeds SEepoch is,000 and 19,000 RPM. 

2. There is a region of instability in the RPM control 
between 17000 - 18500 RPM which is caused by either 
instability in the flow of water through the waterbrake or 
by the inability of the controller to control the turbine 
Output characteristic in this range. 

3. Operation of the rig is satisfactory outside the 
range of the instability. 


4. Variations in Pi , With stage pressure ratio held 
O 


constant, have no measurable influence on the interstage 
or stage performance parameters. 


5. Variations in Ty , With stage pressure ratio held 
° 


constant had no measurable effect on stage performance. 
Changes experienced in stator nozzle pressure ratio and 
stator exit pressure ratios (and consequently, loss coeffi- 
Clients) were directly related to condensation effects 
beginning in the stator nozzle. 

6. Approximately 30 scans of stator and rotor torque 
and axial and closure force measurements were required to 


produce a repeatable average for a given data point. 
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7. The pressure drop across the flow nozzle, are 
has the greatest influence of the three input measurements 
on the calculated value of Weef’ 

8. To ensure accurate computation of RAF, a careful, 
accurate calibration of the force measurements must be 
conducted before each run. 

9. Pressure ports 11 (hood pressure) and 21 (stator 
exit hub pressure) exert a disproportionate influence on 
the uncertainty in calculating the resultant axial force 
(RAF), stator exit pressure, and stator and rotor loss 
coefficients. An increase in accuracy in measuring these 
two pressures would result in a significant decrease in 
the variations found in each of the above-mentioned calcu- 
lated parameters. 

10. The stator exit hub pressure on Port 21, is 
unsteady and exerts a stronger influence on the uncertainty 
in RAF than does the tip pressure on port 24. 

ll. A pneumatic damper on ports 21 and 24 would do 
more tO improve the accuracy of the calculated values of 
the resultant axial force, interstage pressure, stator and 
rotor loss coefficients than any other change which could 
be made. 


12. A nearly constant uncertainty interval exists in 


So @re=.0300 and in oR Oe Uo tO... 08. 
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VI. RECOMMENDATIONS 


l. The speed range of 17000 - 18500 RPM be avoided 
until the waterbrake instability in this area is corrected. 
2. Further tests be conducted to attempt to improve 
the waterbrake stability by perturbing the internal 
geometry. 
3. The pressure tubing inside the hood be replaced 
with high temperature plastic tubing to allow the use of 
temperatures up to 710 °R. 
4, A minimum total temperature of 680° R or higher 
be adopted to avoid condensation effects when interstage - 
data is taken. This can be done by draining the aftercooler. 
5. A pneumatic damper be incorporated in the port 2l 
and port 24 pressure lines to improve accuracy in calculating 
interstage parameters. An improvement in speed stability 
will alleviate, but not eliminate this requirement. 
6. The investigator should check the calibration of 
the AP oz transducer before and after each data point to 


@maeure accuracy in er 
7. The investigator should cooperate with the test 
engineer to verify the calibration of key measurements by 


following, at a minimum, the steps listed in Appendix E. 
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8. Merging the acquisition and reduction programs 
should be considered and implemented if found to be 


peacticable. 
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TABLE I 


TESTS CONDUCTED 


TURBINE PERFORMANCE TESTS 


Run Number TE, Cn) 
8 639 
9 643 
10 635 
ot 639 
1 574 
13 574 
14 602 
LS 676 


RPM SAMPLE TESTS 


Run Number 


Pty 


Pressure Ratio 


("Hg) 


SG 
oo 
oe 
74.5 
Sec 
Boo 
60.0 


Sys) as, 
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Pressure 
Ratio 


RPM 


10040-18870 
TOO -~ 8750 
11077-19346 
toe oO Oa? 
2214-19046 
PAviGo~ 20009 
Oe ao 7 


12081-19940 


RPM 


10600 
S00 
IIL SO 
11600 
11600 
13000 
14100 


Todo 





9 3.0 16100 


L18, S.2 17400 
til. Baw 19000 
2 25 20000 


*RPM samples were taken on most occasions between 
performance test points. No attempt has been made to 
equate these run numbers with those listed as turbine 
performance tests. RPM sample tests are sequentially 
numbered and correspond to the order of graphs of results 
presented in the report. 


DATA SAMPLE TEST 


One test was run to determine the number of samples 
needed to produce a constant average in the forces (two 
measurements), and torques (two measurements). The test 
was run at 11500 RPM, pressure ratio of 3.0, and 100 


samples were taken of each parameter of interest. 
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TABLE II 


STORAGE LOCATION OF RAW AND REDUCED DATA 


TURBINE PERFORMANCE TESTS 


Run Number Record Number File Name 
1 1-10 
2 11-20 
3 erased 
4 21-30 Raw Data in file 
5 31-40 "RAWDAT" 
6 41-50 
7 51-59 Reduced Data in 
8 erased file “REDDAT" 
9 erased 
10 60-68 
11 69=76 
PZ 77-82 
iS 83-90 
14 91-98 
15 99-106 


Pele SAMPLE TESTS 


# Slingers RPM Record Number File Name 
i 11500 1-3 
1 11500 4-6 
it 11500 7-9 
gu S00 10-12 Datayseored in 
1 11500 13-15 file "RPMSTO" 
2 11500 16-18 
2 11500 19-21 
2 10400 22-24 
2 13100 25-27 
2 19000 28-30 
Z 13100 31-33 
Z 14100 34-36 
Z 15100 37-39 
Zz. 19000 40-42 
2 17100 43-45 erased 
2 19000 43-45 
2 20000 46-48 
2 19000 49-51 


UNCERTAINTY MEASUREMENTS 


Run Number Record Number File Name 
10 1-9 
el 10-17 
eZ 18-23 Stored in "DATAZ", 
is 24-31 for performance 
14 32-39 test data. 
15 40-47 
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14330 


Lilo 


15060 © 


15846 
15928 
15046 
15829 
15106 
1SusuV 
1Steb 
15200 
15120 
1Siage 
15146 
151606 
15169 
15246 
15260 
15200 
15140 
15148 
15048 
14368 
15009 
15186 
15198 
151060 
15464 
15080 
15106 
15030 
15646 
151640 
1St2¢ 
15168 
15136 
15126 
15126 
15146 
151368 
15168 
15109 
15126 
150690 
15v6y9 


RPM SAMPLES AT 15,000 RPM 
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TABLE XV 


DEFINITION OF COLUMN HEADINGS IN TABLES XVI TO XXXVII 


SYMBOL 


e=DELT 


DLE 


vl 


V2 


VAL 


VA2 


Al 


A2 


Bl 


B2 


DEFINITION 


Stage AT calculated from the 
measured HP. 


(DEL T) = (C=-DELT) 


Velocity at Stator Exit Plane 
(ft/sec) 


Velocity at Rotor Exit Plane 
(ft/sec) 


Axial Velocity at Stator Exit Plane 
(ft/sec) 


Axial Velocity at Rotor Exit Plane 
(ft/sec) 


Tangential Velocity at Stator Exit 
Plane (ft/sec) 


Tangential Velocity at Rotor Exit 
Plane (ft/sec) 


Mach Number at Stator Exit Plane 


Axial Mach Number at Stator Exit 
Plane 


Mach Number at Rotor Exit Plane 


AzbalsMach Number at Rotor Exit 
Plane 


Flow Angle at Stator Exit Plane 
(degrees) 


Flow Angle at Rotor Exit Plane 
(degrees) 


Relative Flow Angle at Rotor Inlet 
Plane (degrees) 


Relative Flow Angle at Rotor Exit 
Plane (degrees) 
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SYMBOL 

Boor ZL 

ZSTH 

ZR or 23 

ZRTH 

ZR* 

ZI 

a 

esa 

SiIPR 

m.P. 

RTM or DYNA Q 
STM or STATOR Q 


AXF or AX FORCE 


eer Or CL FORCE 


MW-DOT 
EzO 
PO 
BHD 


eo LP 


1 


=-HUB 


P-TIP/PTO 


TABLE <ve cont’ d) 


DEFINITION 

Stator Loss Coefficient 
Theoretical Stator Loss Coefficient 
Rotor Loss Coefficient 

Theoretical Rotor Loss Coefficient 
Rotor Carry Over Loss Coefficient 
Rotor Incidence Loss Coefficient 
After Expansion Loss Coefficient 
Pressure Ratio Pt/Po 

Stator Pressure Ratio 

Horsepower 

Rotor Torque (Moment) (in-Lbf) 
Stator Torque (Moment) (in-lLbf) 


Stator Axial Force - Force 
Capsule (l1bf) 


Closure Plate Force - Force 
Capsule (l1bf) 


Computed Mass Flow Rate (lbm/sec) 
Total Pressure (in Hg) 

Total Temperature (°R) 

Hood Pressure (in Hg) 


Static Pressure at Stator Tap 
Toes Seaeitas Ho) 


Computed Stator Exit Pressure 
(in Hg) 


Static Pressure at Stator Tap 
Baoe #3. (in Ho) 


Pressure Ratio of Pip’? to 
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SYMBOL 
P1/PTO 
P~HUB/PTO 
KIS 


TURB RE 


BELT or DEL T 


BTA 
DELTA 


THETA 


RHP 


RMW-DOT 


RRTM 


RSTM 


RN 
RTH 


Beer 


DEFINITION 

Pressure Ratio of Pi/Pr 
Pressure Ratio of Ean Eos 
Isentropic Head Coefficient 


Reynold's Number at the Stator 
EnErance 


Temperature Difference Across the 
Stage = Tt, = aoa 


Total to Static Efficiency 
Pto 
Pressure Ratio,é6= 5 





(PB 


oe ref 77° 22) 


Temperature Ratio, 9= rT 
(T of — 518. 7°R) 


Referred Horsepower 


Referred Computed Mass Flow Rate 
(lbm/sec) 


Referred Rotor Torque (Moment) 
(in-1bf£) 


Referred Stator Torque (Moment) 
(in-lLb£) 


Referred RPM (RPM) 
Theoretical Degree of Reaction 
Effective Degree of Reaction 


Resultant Axial Force (Computed) 
(Be) 


Rotor Speed (RPM) 


Pressure Ratios Corresponding to 

Tap Locations in Figure 6. Pressures 
are referred to Stator Inlet Total 
Pressure (Py) 
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APPENDIX A 


COMPUTER PROGRAMS 


A.1. ACQUISITION, REDUCTION, AND BATCH PROCESSING 
PROGRAMS 


The acquisition program used to acquire data in 
Ref. 1 was entitled "SCATTR". The name was changed to 
"SCAMOD" and the changes to the program are listed in 
Table A-l. Listings of the current acquisition and reduc- 
tion programs are given in Tables A-II through A-Ix. 
Programs TTR7, TTR8, and TTR11 are currently used to batch 
process raw data, reduced data, and uncertainty data, 
respectively. Listings of these programs are contained in 
Tables A-X through A-XII. Reduction programs were also 
modified, and the changes are listed in Table A-I. A new 
program, "UNCERT", was added to the reduction sequence to 
calculate uncertainty in the loss coefficients, RAF, and 


P1/P The formulas used in program "UNCERT" were 


t- 
BP eained from the uncertainty development in Appendix C. 
A current listing of the reduction sequence is presented 
in Figure A-l. It will be noticed that the smoothing 
procedure using TTR1 and TTRI has been deleted from the 
reduction sequence. (Examination of Figures 31 and 32 
clearly showed that addition of the "slingers" to the 
waterbrake assembly resulted in much less noise in the 
calculated values of Pie yielding a smoother variation of 
P,/P 


than was seen in the data gathered in Ref. 1.) 


La) 


ee 





Line 590 in TTR2, which allowed setting Woe = 
1.02/lbm/sec = constant, was deleted since it had never 
been shown that the value should be constant at all test 
conditions. It will be noted that except for run 12, the 


W for runs 10-15 was fairly constant at around 


ref 
1.025 lbm/sec +/-.005. 

The storage and batch processing programs (TTR6 and 
TTR8) were modified to print out engineering data and 
other parameters. Some consideration should be given in 
the future to expanding the storage array to include 
engineering and other data. There is room for more data 
elements per record.than is presently stored. For example, 
the Labyrinth flow rate, flow nozzle pressure, temperature, 
and pressure differential should be stored. 

It was necessary to "scratch" the acquisition program 
at the end of each data point in order to reduce the data 
Just taken. It is recommended that the feasibility of 
merging the acquisition program with the reduction sequence 
be studied and undertaken if possible. The only limiting 
factor is the number of variables left from the reduction 
sequence which can be used in the acquisition program. A 
current list of the variables used in the reduction sequence 
@mey is listed in Figure A-2. 

Table A-XITI contains a list of all new variables added 
eemene reduction programs. A list of all plotting routines 
located on the TTR disc and variables plotted are listed 


in Table A-xIV. 


a3 





The programs used to acquire the RPM samples, VREMSUR, 
and to compare calculated and measured temperature drops, 
"TERRY7", are listed in Tables A-XV and A-XVI respectively. A 
list of all programs used for acquiring, modifying, or 
reducing the data presented in this report is given in 
Table A-XVII. 

On runs 13 and 15, variations in the recorded voltage 
values for a pressure port measurement indicated the 
pressure line had developed a leak. Accordingly, the 
reduction program was altered to compensate for this 
Situation. Table A-XX lists the pertinent runs and points 
applicable, pressure port, and action taken in the 
reduction sequence. The voltage values are listed in 


Tables XX and XXII for Runs 13 and 15, respectively. 


A.2 FOURIER ANALYSIS PROGRAM 

The program used for Fourier analyzing the RPM data 
was program "AUTOST" which is contained on tape and asso- 
Ciated with the HP9845A computer system. The "DATA INPUT- 
TIME DOMAIN DATA" subprogram option was used and had to 
be modified as listed in Table A-XVIII. The data files 
containing the RPM data analyzed are recorded on the 
"TPL Library #1" tape in the NPS Turbomachinery Laboratory. 


The data file names are given in Table A-XVIII. 
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PROGRAM 


SCAMOD 
(see note 1) 


TABLE A-I 


LINE ACTION 

100 Added array E(2) 
L0 added 

160 added 

180 added 

190 added 


210/240 added 


350 and 
1710-2280 
(see note 2) 


860-990 delete 


temperature 


averaging 


870-1690 Recorded 
channel 
numbers 
changed 


(see Appendix B) 


le as 


CHANGES IN THE ACQUISITION AND 
STORAGE PROGRAMS 


COMMENTS 


Averages polnt. 
24 pressure 


Enters run #, 
day, month, year 


Takes multiple 
DvoOints without 
restarting program 


Loops program for 
multiple points 


assigns point 
number for data 
run 


AsSigns point number 
when auto-storage 
not selected 


Records 46 pressures 
On scanivalve 
system 


extraneous 


Changed channel 
setup 





rR IB 


ER 2 


PER 3 


UNCERT 


TABLE A-1 (Cont'd) 


1050 and 
22 20— 29 60 
50 


40-50 


443 


451 


632-634 


a 56 


770-1080 


POU 20 


1060 


PO F07 aE 0), 


PrCO7 foo; 


1190 


iO 


590 


eee7 


580 


All 


Port 24 
averaged 


added N(20) 


added 
Councer 


nod2i£1ed 


modified 


added 


added 


modified 


modified 


changed 


modified 


changed 


deleted 


added 


modified 


added 
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Averaging 
sequence 


“UNCERT" program 
St@Gbage array 


Record # for 
storage of 
"UNCERT" data 


channel number 
assignments 


Channel number 
assignments 


turbine outlet 
ECnPperatuLe 


pressure conversion 
EO. Ins. Aa 


scaling factors 
assignment 


scaling factors 
assignment 


more stable port 


channel 
reassignment 


averaged channel 
used 


eliminates 
constant W 
ref 


facilitates 
multiple single 
point reduction 
without smoothing 


Ghanins sto. UNCERE” 
for uncertainty 
calculations 


calculates 
uncertainties 





ick 


PRO 


0k / 


TTR8 


er 1 


TABLE A~1 (Cont ‘d) 


60 modified channels 
reassignment 
85 added prints T 
O 
Sea) added prints calculated 
WwW, W 
ref 
105 added turbine stator nozzle 
Reynolds Number 
DELnecout 
534-538 added Prip/ Pe ; Pi/P, ; 
O O 
and Paup’ Pe printout 
470-510 modified different channels 
used 
520-562 added stores 19 other data 
parameters (see 
note 3) 
690-692 added Chains tom@inrip for 


multiple point 
reduction with no 


smoothing 
270-330, modified channels reassigned 
480-780 
348 added program loop more 
efficient 
410 added "Table" format 


Completely renumbered and reformatted to 
conform with 19 additional parameters stored 
aha ye Or 


All added Batch processes the 
“UNCERT" data 


ae 


TABLE A-1l (Cont'd) 


Note 1: Due to the extensive nature of the modifications 
to "SCATTR" (see Ref. 1), the program was changed to 
"SCAMOD" and renumbered. Major additions to the program 


are listed in the table. 


Note 2: These lines replaced lines 110-210 of "SCATTR" 
program (see Ref. 1) due to replacement of the paper tape 


system with on-line data acquisition. 


Note 3: The nineteen additional stored items, beginning 


with line 520, are: 


Line Program Variable Parameter 
520 Z9 W 
BAG v9 W 

ref 
530 Q8 Re No. (Stator) 
OZ Hoe Horsepower 
534 M5 Rotor Torque 
536 M6 Stator Torque 
538 F4 Stator Axial Force 
540 F5 Closure Plate Force 
542 R2 Resultant Force 

on stator 

544 On Py 

Oo 
54 

6 ve Te 

548 Os Hood Pressure 
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550 


2 


554 


26 
poo 
560 


6 2 


TASLE eA be (cont d) 
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Sf 


R8 


04701 
P2 
R8/Q1 


2 


129 


Static 
stator 


Stator 


Static 


pressure at 
tape iip. tor 


Exit Pressure 


pressure at 


stator tap "HUB #3" 
Porp/ Pt 

O 
P,/P, 

oO 

Puup/ Pt 

O 
Temperature difference 
across the stage 
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TABLE A-XVII 


PROGRAMS USED FOR ACQUIRING, MODIFYING, 
OR REDUCING DATA 


PROGRAM USE 


SCAMOD Acquires raw data for turbine 


performance tests 


oO RI B Recalls raw data from storage. 


Assigns scaling factors. 


ie RZ Calculates W W, W and 


LAB. ref’ 


calculates parameters in stator 


control volume (So, etc.) 


oR 3 Calculates parameters in rotor 


control volume (he etc.) 


UNCERT Calculates uncertainties in 
RAF, Pi/Pogs Cor and op 
TTR4 Calculates stator velocities 


and theoretical loss coefficients 


mrR5 Prints out results for each 


reduced data point 


TTR6 Stores reduced data 


Zat 





eR 7 


TTR8 


rR 1 


RPMSUR 


BERRY 7 


Tabulates 


Tabulates 


Tabulates 


Takes 360 


raw data 


reduced data 


uncertainty data 


RPM samples ata 


rate of 330/min and stores 


them on a 


set of 3 sequential 


record numbers 


Compares actual and computed 


temperature differences across 


the stage based on horsepower 


calculations 
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GABLE A=XVIII 


WAVEFORM ANALYSIS PROGRAM MODIFICATIONS 
TO ACCEPT TINE DOMAIN DATA 


Statement # 


2520 (Current Program Statement OK) 


2521 (Begin Modification) 


Ee 30 


2531 


B32 


2533 


2534 


2540 (This Program Statement OK) 
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Statement 


OP Ee art k 


POR] i= =) TO 12S 


READ #1; Real (TI) 


NEST 


Powel = 7 TO 128 


READ #1; Imag(IL) 


NEXT I 


temeeony o1.,Ll)="N" 
THEN REWIND BS 





TABLE -A=ALX 


DATA FILES AND PROGRAMS USED FOR 
FOURIER ANALYSIS 


Data/Program 


15000 RPM (DATA) 


19000 RPM (DATA) 


40 Point Power Corre- 
lation Data 


FIXER (Program to take 
out the 'DC' components 
from raw RPM data and 
store on a different 
data file) 


Note l: 


Where Stored (Note 1) 


oTERRY 2°’ (Raw Form 
256 Data 
"AVG2" (Averaged 
Z56 Data 
See Note 
"TERRY1" (Raw Form 
256 Data 
"AVG1" (Averaged 
256 Data 
See Note 
"TERRY 3" 


TPL Library Tape #1 


Points) 


Form - 
Points) 
Z 


Points) 


Form - 
Points) 
2 


(Before referring) 


All data files listed in this table are stored 


on TPL Library #1 Tape (Not to be confused with 
the H.P. Library tapes.) 


Note 2: 


The data could only be analyzed in powers of 2; 


therefore all 359 points could not be used and 
Only the first 256 points were analyzed on each 


RPM data set. 
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Run 


3 
it 3 


ES 


Note l: 


TABLE A-XxX 


RAW DATA ANOMALIES 


Point Pressure Port 
ui 8 
8 8 

2-8 aE ee 


(see Note 1) 


AGELon 


Deleted A(8) 
from line 1080 
ier RIB for 
points 1 and 8 


Used point 1 value 
fom all points in 
line 1090 of TTRI1B 


AP is recorded on channel 27 through scanner 


OZ 


number 1 (see Appendix B). 


Zou 





SCAMOD 





UNCERT 


TTR4 


ik 
TORO 
7 
a” 
a 


a“ 


a oa” 
7 


7 
a 
Ta aay = Aer: 


FIGURE A-l. Acquisition, Reduction and Processing 
Sequence 
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FIGURE A-~2. ARRAYS AND VARIABLES USED 
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APPENDIX 5 


TEST RIG MEASUREMENTS AND 
SCALING FACTORS 


The test rig measurements were obtained uSing the test 
instrumentation of Figure 5. The 46 recorded scanivalve 
ports and corresponding quantities are given in Table B-I. 
The channel numbers and corresponding quantities for the 
15 other recorded parameters are also listed. The scaling 
factors applied to the raw data voltages to obtain the 


proper engineering units are given in Table B-II. 
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TAREE 2sa1 


TEST RIG MEASUREMENTS 


Scanivalve Port Number 


(Scanner 1) 


Wo WA UW & WD - 


29 through 46 


Pressures 


Quantity 


Atmospheric 

Calibration Reference 

Flow Nozzle (upstream flange) 
Flow Nozzle (downstream flange) 
taeyrinth Plenum (Fig. 2) 
Stave Potal Inlet (Fig. 3) 


Hood 

Seaconrulanegi (Fig: 7) 
Stator Tap #2 (Fig. 7) 
Secor lap to (Fig. 7) 
Seatoretap #e"(Fig. 7 
Seaton lap +o (Pig. 7) 
Stator Tap #6 (Fig. 7 
Stator Tap #7 (Fig. 7) 
Stator Tap #8 (Fig. 7) 
Stator Tap #9 (Fig. 7) 
SEcOumeLee nub (Fig. 3) 
Stator Tip (above Tap #8) 
Stator Tip (above Tap #9) 
Seateer Gxtearrp (Fig. 3) 
(See Note 1) 

Atmospheric (Spare) 
Calibration Reference (spare) 
Labyrinth (spare) 

Stator Inlet Static 
Rotor Shroud (spaced at 
l inch between centers along 
shroud inclined face 
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Scanner Scanner 
Number Channel 


2 00 
Ol 
O02 
03 
04 
05 


06 


TABLE B-I (Cont'd) 


Other Channels 


Quantity 


Air supply temperature 

Flow nozzle temperature 

Inlet plenum temperature 

Closure volume temperature 

Stator total inlet temperature 

Turbine outlet total 
temperature 

Labyrinth temperature 


RPM 

Stator axial force 
Closure plate force 
Dynamometer torque 
Seacor Eonque 
Atmospheric pressure 
Flow nozzle pressure 
Flow nozzle AP 


Note 1: This channel is sampled 30 times and averaged 
prior to storage to reduce unsteadiness. 


BENE 





Parameter 


Pressure 
(all scanivalve 
ports) 


Temperatures 


Axial Force 


Closure Force 


Dynamometer 
ZOorque 


SeacOr Torque 


Atmospheric 
Pressure 


Flow Nozzle 
Pressure 


Flow Nozzle AP 


PaABbiles =i 1 


SCALING FACTORS 


i 
(Scanner 


00-06 


(Scanner 


16 
(Scanner 


20 
(Scanner 


21 
(Scanner 


23 
(Scanner 


24 
(Scanner 


25 
(Scanner 


26 
(Scanner 


24 
(Scanner 


Channel Number 


#1) 


# 2 ) 


#1) 


#1) 


#1) 


#1) 


#1) 


#1) 


#1) 


#1) 


Zou 


Scaling Factor 


Conversion to 
IneHg x LOO0 


Polynomial regression 
FOEmuUra tO °F e+ 460 


EOgrok 





APPENDIX C 


UNCERTAINTY ANALYSIS 


An uncertainty analysis was conducted on four key inter- 
stage parameters, 1.e., stator exit pressure ratio (P,/Pt.), 
resultant axial force (RAF), stator loss coefficient (2), 
and rotor loss coefficient (cp) - The analysis was done in 
accordance with the principles contained in Ref. (16). The 
following formulas were used repeatedly to develop the 


propagation formulas for each intermediate parameter of 


interest: 
Note: "U" denotes uncertainty interval throughout this 
appendix. 
Up dln es wv 3 d1ln Ww. YW D ay ae 
Desired form: 7 “Sink V,? (Th R wy } 


where ue and w, are independent variables and R = f(v,,w,). 


merce R = x + y then 








en x Xx aa SOE vs 
dln R Xo ain ok nee 
and 
1/2 
im 2 2 
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elfen — lf — yy, then 





U n U 
ee a 
R 1 - yn y 
n 
3 ten = x + yy, chen 
UR = x Uy 2 ny" ae 1/2 
=e = cee ee a y J 
x+y x+y 
4 Mie =axy ++ Z , then 
oR = ((2¥e SEP fa caxy” : Wy) 2 E (mz el 
R R R 4 R 
5 tf-R = V , then 
ed 
R V 
pee Li R = ax + by, then 
iy 2 
Ue eenamO )s beu ) : 
R x 4 


The following development is divided into five sections. 
In each of four sections, one reduced parameter uncertainty 
is calculated. A simplified analysis section follows. 

Table C-I contains the definition of terms used in the 
analysis. Symbols used in this appendix are consistent with 


those used in the reduction program in Appendix A. The 


259 





assumed uncertainty intervals in the raw data parameters of 
interest are given in Table C-II and are consistent with those 
given in Ref. 1. The point chosen for analysis was run 10 
point 1. Pertinent raw data, and intermediate calculation 
data are given in Table C-III. The uncertainty approximation 


formulas developed here are summarized in Table C-IV. 


C.1 RESULTANT AXIAL FORCE (R2) UNCERTAINTY 
Using Eqn C(4) in Reference 1, the resultant axial force 


1S given by 


R2 = F4 + F5 + 0.491159(F3 - F2-F9') 
where 

FO' = FOG + Fl + AB + AD + BS + BE + B7 
and 

FZ = N8(0.25) + N9(0.50) + 01(0.50) 


The uncertainty in R2 is then calculated by first calculating 
the uncertainty in each of the terms in the above expressions. 


These calculations follow: 


1/2 
Minos OSy- + (.5x .04)- + (.5« .04)°] 


G 
| 


FQ 


= .03 


240 





since, 


Fl = 0O.51(83 + S4 + §5) 
TG seamen <4 ue) ~ au * al 
Fl $3 S4 S5 
Baye 
= 0.51[3(.04) *] 
Posse 
and since 

A8 = 0.51(S6 + S7 + S8) 
U Oem Ul se CU a 
Ag Sé o7 $3 


3 eZ 
Oris (.04) ~] 


Ul oS 


Also, since 


AQ = 0.52[(S9 + NM + N11] 
2 
z 2 2 2 
7 2 


DeSee04N-j 


= .036 
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Similarly, 


B5 


G 
it 


EBS 


B6 = 


Bo 


Omeozives a0. o2NGe 4.0. 53N4 


ec 2 2. 1/2 
2 


[(0.52U,, + (0.520, eo) 5S Uae | 


3) 


2 2 2 V2 
mms - eer oes2. (04) ° + (0.53 .04)~*} 


-03'6Z5 


Ore eatoet 4. 50999(N6 + N7) 


72 


g 


2 


2 
Ne? 


GOre muerte. 80999 U to ee OO oe Wee) 


N5 N7 


D > ie 
COs? bx 04) + (4.80999 x°.04)~ «x 2) 


22750 


Note here the order of magnitude difference in the uncertain- 


pepper a rm me rm SS SSS 


mmm em mm i me 8 


Continuing, 


B7 


By, 


Therefore, 


it 


1.897 x 04 


oo) Uo4 = .07588 
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2 Z 


2 


Cc 
tl 


+ (Ugg)? + (Ugg) 


2 2 21/2 


eee (Ul gee (Cae) 2} 


[(.03)° + (.0353)* + (.0353)* + (.036) 2 


2 2 2,12 
+ (.03625)* + (.2736)7 + (.07588)2] 


= 2943 
Samilarly, 
F2 = 47.66118069 x R8 
Una = 47.66118069 Ung = 47.66118069 x .04 
= 1.9064 
Howe  o2.51589456 <x 03 
Un = 82.51589456 Uo3 


= 82.51589456 x .04 


= 3.3006 
finally, 
Wie = nus) ie (Uns) + (.491159 Un3) 
, 2 2 
+ (~.491159 U,,)" + (~-491159 Ung.) 7] 
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Ge emg (451159 x 3.3006) - 


G 
il 


bez 
2 


+ (-.491159 x 1.9064)* + (-.491159 x .2943) 7] 


OGL Oe 


meee re the predominant effects of O03 and R&, i.e., the 


Meeaepressure and stator exit pressure “hub #3" respectively. 


Therefore, the potential % error in the resultant axial 
force is equal to 
Te 1.9966 


R2 x 1003 = 1S0n79 1003 = eos 








C.2 PRESSURE RATIO (P2) UNCERTAINTY 
Since Ql is the average of five pressures (Ports 6- 10, 


see Appendix B), 





D i 2 

Yo1 mi ea UP ch eee =  .275(.04) 

= 20179 an He 
= Ql 
Q2 = 39.93 
Therefore, 
noe) om | 50179 = 0003 
Q Q 59.622 : 
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Similarly, 


From this point on, Wee 





_ m2 N72 
T3 = (Syg77) 
vP3 el Yp2 eee =. 0039 
ES Ze 2 636.40 R : 
wy = Wref Q2 
3 
U Oe U iy 2 
we _ r( ref)2 + ( Q2)2 ae 13) 2, 
wd Woe QO2 Ts 


f (referred flow rate) 1s assumed to 


be constant at 1.02 lbm/sec. Hence 


and 


Bontcinuing, 


Uae ae 





a OO soi- + 000s) <1 = .00049 
U U 

= a 00049. 

Vl = 109.62 SQR(T2) 

U GB 

tet) 
Sr = 5 ES = «-00039 

ee ee 1,2, ,. 91,2,°"° 

v7 Ww val O1 


Za 
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U 
mee = [(.00049) * + (.00039) * + (-.0003) *] 
= .0007 
_ M6 
=) wAnenA gas 
Yves _ (M6) 2 ae Mwa 2,77? 
VG M6 w4 
2 
: 99 ,2 _ 2 
l (sEy-Gp? + (-.00049) 7] 
= .0027864 
_ Vg 
XG = VL 
me fly? (- vn)2,°7" 
XG V2 V1 
2 p hve 
= [(.0027864)~ + (-.00039) “J 
eo Ge 
aa == Saas l(V7)* x (1 Segoe 018 (A+) “] 
Let 
m= 1 = (xg) -- 
pence X89 = .5435 (using data for this point), 
Z U 2 
OM = (00281) =e O26 
1 = (X@) 1 - (.5435) 
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Let A = (V7)°L, then 


[2 te)! 


Ya van mee 
A v7 i 


Why 


[C2 x 0007) * a MOOee oN 


= .00274 
\ R2, 2 
Let B = 018 (AF Then, 
ee? 4 gi,2,’”* 
B R2 ep 
72 


2[(.0153)* + (-.0003) 7] 





= .0305 
mnerefore, 
aoe = .Llov(A— B) 
“28 p (2 tO 7A) aye fe eens a 
ZY A- B A a. B 3B 


From the definitions of A and B, 
A = ca and B = 018 (— 
: O1 


Ser this point, 
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ve S857 1 
2 

on - == 7045 

Foe] 130-79 

ies NGa2 

Hence, 
Pa ees (7e45) = Joege 
_ ieee. 

Me EMEC gc 42 do 
ZB mee ndc6S 


~ 1353 


meron (09809) 
mols — 20060) 


Wee 


a 


) x .0305) 


Note here again, the predominant effect is from R2 in the 


© mre prem mm rrr ae mS ES SS ES 


Emeertainty error of 2M. Continuing, 


_ R2 
BY = .0558 OL 
Bee a? oo. 912,77 
BY R2 Ql 
1/2 


Taos = 0s) —] 


Now 


P2 = BS + SOR((BS)* - 2B) 
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Oba 3 





Mec = (Bg)* - 29, then 


Uc 2(Bg)* UBg.2 20 V2G, 2 1/2 
“Cia aa Bo? oe Oo aH ] 
(BY) - 29 (BY) - Z¢ 
For this point, 
cf 
Wgmeeeseeie267 < 10 * 
therefore 
U 2 
= - ee 
Grieoane = (5040. 10°") 
aan 1/2 
+ (= —Sdge7 2074 3353992; 
(Opa 205649 1087) 
0321 


—_ rrp mm mm se Se 


0 mm me mm me See Se 


P2 = BS + VC 
ey seal 2 
m2. (Be 2, ZS Mey 2401" 
P BY+SOQOR(C) BY BG +SOR(C) C 
Here, 
a 2 -4 
eee ee 2aye = {5.49% 10 ) = .01445 


Vena) 202. 
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Therefore, 


=(.1202) 1/2 


pee ae 202 


a2 mee (.0321))*] 


: 2 
p> ot OU (qty ez 8 193)) + 


Opel 1 


am mmm mer er rm me eS ee 


See OLATOR LOSS COEFFICIENT (S,) UNCERTAINTY 


The calculation of the uncertainty in stator loss 
coefficient requires first the calculation of the uncertain- 
ties in the intermediate variables defined in the data reduc- 
tion program. In the calculation of the axial component of 


Ene velocity, 


P P2 

x1 = 3.487562189 77 
pate P22 Uy7 2 +72 
Se an (55 a = TF? ] 

= 2 2 = 

= Owe e+ (=.0007)")-- = £0111 

Then 
— t t 2 Z 

Xl = = X1l' + SQOR((X1')~ = (XM)~ + 1) 


MeeeD = (X1')* - (x9)* +1. Then, 
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> prea) — een 2 2 (xg) 7 Ug 21/2 
— ~ i —_* > 5 (a7) ) 1 or oS 5 (xq)? ] 
(X1')° — (XG)~ +1 (X1')~ - (Xf)"~ +1 
fer this point 
a 2407 
Del: = 3.487 (Seay) = 2.37 
X89 = .5436 
Therefore, 
U 2 
ws = (542 ___(.0111))° 
(237) (. 5436.) + 2 
cao Ife 


+ (- — 52 '3) (00281) ) 1 
Crs (2.5456) + 1 


SOLO? 2 


e 


The predominant effect here is from P2 and therefore from R2. 


mance Xl = SOR(D) - Xl', 





eee 
a. 7 o “D2 eae X1' i 
5 awe, D s/o Xl 
mee this point, 
Zz Z 
Dee) (2237) = (2.5436) + 1 = 6.323 


so that 


Zou 


U 106.323) 1/2 
x1 eee 61972) ) - 
x1 173 (.0 
icpeconeo a (2537) 
1/2 
C= oath enone + | 
(euseeNe = (2.37) 


a. 206 


mean, the prime contributors to the uncertainty in the axial 


Note also the high errors introduced through taking differ- 
ences between calculated quantities as in the denominators 


of the bracketed terms in the Xl uncertainty formula. 


Bentinuing, 
K2 = SQR((x#)* + (x1)?) 
et 
U 2 U 
Meixe)?+ (x1)? + 2 - (2x) _ ae 
(XO) “ + (X1) 
ema? 
Pimmoueme > xi | 
(Xf) ~ + (X1)~ 
Ber this point, 
Xl = .144224719 


Therefore, 


25:2 


2(.5436) 2 
(.5436)° + (.1442) 


tl ey 


= [{ =(.00281))* 


2(.1442)2 whe 


> ( 5 
(.5436)° + (.1442) 


2 
5(.2506)) 7] 


eee 03338 
Mee = SOK E 
"x2 — rE = 0167 
X2 Carl ; 


Mempencough its effect on P2. Finally, 


2 
21 = t. = 1 (<2) 55 
1 - (p2)° 
Met 
ems = p27-~’ 
Then 
EG) . ee, 22 
G e287 2 
287 1/2 
_ pee 2287022427) ae ne 
1 - (.2427)° 
=) enyeee 


255 





Let H = a Then 
CH = ne? Oe Z SS 2)" 
H x2 G 


Wy 2 


meee olla?) ~ 4+ (2) 00leed) —) 


303395 


mmr mm mre eg mee | em 


m@meee2, the overriding importance of the uncertainty in R2. 


montinuing, 





Zi = 1-48 
2 
zy 1-H # 
Pea (x2) 
G eee 
for this point 
X2 = .562373442 
Therefore, 
2 
a= C502) ~-- = .9478 
iene 2427) - 


and 
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—ee eS Se mmr Se Sm ee 


Zemene Major contributor to the uncertainty in the calculated 
Maeerstage pressure and consequently to the stator loss 
coefficient. The uncertainty in R2, however, is in turn 
determined primarily by the uncertainties in two key pressure 
measurements. The measurements corresponding to computer 
program parameters Q3 and R8 are those on ports ll and 21 
which are connected EOmENe IiI@eGe ana cto the closure plate 


volume. 
C.4. ROTOR LOSS COEFFICIENT (Zp) UNCERTAINTY 


The uncertainty in the rotor torque is obtained using the 


measurement uncertainty: 


U | 

Pope o7an-lbe 

mer 446.90 = 0.000955 
Then, 

Kaue=ge4 xg — 1235 — “to __ 

Wa (VL) 

Let 

; = _sMS 

wW4 (V1) 
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U U U U / 2 
— M5, 2 _ ~W4, 2 _ Be 
a [ (sae) + ( wa + ( v1 ] 
2 2 21/2 
= [(0.00095)“ + (-.00049)“ + (-.00039) “] 
= ,00114 
Therefore, 
x4 = 5 (__- 98 4X9 pe eee 351 pee; 
aC ~984X8 - .2351* XG “984xX9 —= 52351 ° I 
_ M5 o 446.9 2 
ED - Watvl) ~ (Clos7a)(o7es.3e))60OU C7 BLS 
Let 
J = U3(X@) - U4(X4) 
and 
Kee = U3(X?2) 
L = U4(X4) 
U U U u/2 
on W232 XG. 2 
—— | L(+) a (<r) ] 
2 21/2 
= ((,0853) + °(.00281)—] = .0161 
L = (ud “ + (x4 ae 
Cc U4 X4 
2 2.1/2 
= [(.0158)° + (.00696) ~~] = ol 


1 
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Therefore, 


2k 2 1 
giving 
"Se (eK KP 4 (oR 2)” 
a (aman K-L TL 
U3(X®) _ .0365(11077)(.5435) _ 
K ni) TECiCSmSGl5S5oe se Uieee 
- Sos 7Ilia77) (211272) _ 
L = wv4(x4) es 0189 
cen (5506) = -235(2.815) = - .1272 
U 1/2 
i 07949 2 ~(-.0189) 2 
Meee 79490 -(-.0189) 9181!) + (-aagqg = 7[—loTaa) f° 0172) 
= 01338 
Let 
: 2 2 
Pei we = 1d =| (=.1272)- = .9638 
therefore, 
Um — 2(x4)* Uxe 172) - 
Se SE = EEL (006986) 
ene) Lo (S129) 
= .000229 
72a ee =a ) 


21 


CG 


X4 


Z2 


U3 


Looe ee 2429 —=—9(=.0189) = .0984 


i — ee ee “s,2,"” 
- 2J M- 20 J 





9838 2 
l (3933 - 27.0984) (- 000229) 


=2(.0984) ie 


2 
9838 -2(.0984) §°91338)) 7] 


ae 


00230 


mieel. 24355) 


2 
l994(.5435) — 235 (2815) + 00281)) 


Bea ene ES) , ete 
= 


i mORTEaESa SiN aSKay ors) 


.00696 


(l - (x4)2) - 2(U3) (XM) - (U4) (X4) 





Ul 
Wel 
-J3651ND 
U 
ss 
yr mites a OLS 79 
U U iy AF eZ 
Ul, 2 v1, 2 7 oo 2 
l(a + (- VL ] =o 5 79)" + (=,00039) ] 
2O1 56 
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Now, 


U4 


where 


UZ 


Therefore, 


G 


U2 


and 


U4 


eo 


S9 


So 


S9 


OS: 


= .0371(N) 


Ooo 


| 
I 


Ul 


Oa SiG 


+ (225.93 Uo3 Fe 
SOFO3) OF 


0.5(s9) Ys9.2 
S9 +03 89 


0.5(.04) ae 


5 fF 0.5(.04) 
34.2929 


2 
34.2909 | 


»o eG 


= .0008247 


= 17.0582 


04 


17.0582 =~ 90234 


=e? 2547 


Zou 





es 04 


OS «(OTF aaa7 = 00232 
- Pg 
Bl = 38.78(=77e7y) 
and 
Bo mee 17. 14645 
therefore, 
ly BS? 2, - Ovd,2,°"° 
oe P9 w4 V1 
2 2 21/2 
= [(.0008247)* + (-.00049)* + (-.00039) °] 
ZO ones 
X5 = -Bl + SOR((B1)* + 22) 
Let 
7 2 
N = (Bl)~ + 22 


and since 


Bl 4.188649456 


Z2 


- 78699067 
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U a2 


N 2 Bie 22 Pe 
(Bl) " + 22 (Bl) “+ 22 
2(4.188) 7 : 

a  ! 

(4.188) ° + (.78699) 
1/2 
+ (382? _____( 00336) ) “1 
(4.188) + (.78699) 
= .00198 
Therefore, 
X5 = -Bl + VN 
and 
N = 18.33177494 
CxS = (22 SOR(N) oN, 2 aca =p il yore ae 
eS “Bl + SQR(N) N =BL+SQR(N) Bl 
72 
= [ (2-2 08.33) (00198) ) ° 
=4.188+ (18.33) 
is 1/2 
-4.188+ (18.33) 

= .0653122238 


Note that it is the differences in the denominators of each 


mt gr rm mm mm 


crm ee ee SS 


large resultant uncertainty. Continuing, 


26.1. 





Lg 
Let 

O 
Here, 

U4 
Therefore, 

no 

0 

LY 
and here 

X5 

O 


Therefore, 


(x5)* + (x4 - U4)? 


x4 - U4 
14856054 
y Se 1/2 
x4 Yx4.2 -u4 4.2 


l(ya=ua xa? + ‘xa-ga wa? ! 


=o 3 2 
[ (7593 — qa 9G (- 00696) ) 


mee 2 1/2 


ie miei 922!) I 


.009096 


2 


(X5)7 + (0)? 


aUlezgt2751 


ecole O42 
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Bet 


net 


Spey 2, 2lo)* 0, 2,7 
LY (x5) 2 + (0) * i Coo . 
2(.0929) * 2 
= [( 5 x (-0653)) 
(.0929) 7 + (-.2759) 
_ 2 1/2 
+ ( 2 ( 2. = (00909) ) “] 
(.0929)* + (-.2759) 
= .021073843 

Ll = (1 - (x2)*) (= (GS) 

POR Se ogius)) 4 (ud) ~ 
P = i- (x2)° 
U 2 U 2 
vB rier) ~X2 _ 2(.562373442) (0167) 

ea) Meneses) 

EOE doors 
pg 
2 = ~p2(Q1) 
U U U U WZ 
fe ee ee,2.. 22,2, “9,2 
6 (S57) + (57 + (ST? } 
1/2 


=O OGe 2 ce 0111)° + (.0003)~}] 


oe 5 OIG IEE 


os 


meee that the P2 term predominates here again. Let 


too 
. 17.14645 
= Ba(Gl) ~ .242704755) (59.62218) 
= 1.1849 
Pee feddarce Yo _ ee ota 
: Pons © erecie sao ° 
= .0670486725 
Let 
S = PxR 
Then 
U U cea 
= (aes 
S P R 
2 pee 
= [(.0154)7 + (.067048) 7] 
= .0688 
2 C 
T = (x2)? 
U U 
T _ X20 Z 
ew = 2 SS = 2(.0167) = .03339 
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inet 


met 





feel 


2X9 (U3) 


U U i772 
XO, 2 Warez 
pone uae 


D 5 ey 2 
GsOlzie4)) wa. 0158) =] 


-01605 


U 
2—- = 2(.0158) = .0316 


Saree Vy a 


Pemex) (1 = (Oo)? 

Meee peoyereazy i — (21-1849) °°°°’) 
-.034 

(x2)7 = .316 

2(Xf) (U3) = 2(.543565192) (.146240442) 
15898 

(u4)° = (.1485)* = .02207 
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Therefore, 


ml 


i 


Therefore X2 


Binally, 
Z3 
Let 
X 
Then 
Ux 
xX 


ed 52643558 


U Un 


i 2 
(TO) 


-.034 


2 
[(—yqex (0688) ) et 


= lio 


2 
3453 e2Om OU >) 


aoa 


.07664 


U 


U 


se 
“= = 0.0726. 


YT is the most important term. 


7 _ Lg 
eee | orl 
Lg 
Ll 
U UL, 4 1/2 
(apy? + (= SB * 


[(.021)* + (-.07664)°] 
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L/ 


2 


-07948 


“V _V,2 WLW) 2,1/2 
nee | 7 ‘ni we 
36 2 
F753 §-93339)) 
iy 
~-02207 2 
+ (37597 (-0316)) ] 





23 = t = 1- xX 
where 
ee ee 
Therefore, 
8 7 ane “x = (FSS gy) (.07948) 
= «liz 


C.5 SIMPLIFIED ANALYSIS 

Given the predominance of the R2 term in both loss 
coefficient calculations, we can simplify the uncertainty 
calculation at any point by dropping all uncertainty terms 
which have a magnitude which is less than one fifth of the 


R2 uncertainty. With this approximation, 


CG 
ll 


1.9966 lbf (exact for all values of R2) 


R2 
Or 
U 
ees) ss 
RD RD or a given R, 


and the uncertaintiies in Ql, Q2, T3, WO, W4, Vl, V7, V@ and 
X@ will be considered to be negligible. Under the same 


approximation, 


Zoe 





G 
G 


"Bg _ “RZ 
BG ~—_—sC R22 
and 
Yo _ _2(pp)*_ Tra 
Ee (BG)? - zg R2 
where 
C = (Bf)* - 2p 
and finally 
"P2 _  5;_Ba eo? geese *  Ye.2,* 
P2 B@+SOR(C) BP B@+SQR(C) C 


where again BY and ZY must be obtained from the computer 


reduction for a given point. Then 


als) | | 2 
a" P2 
Sy 2 ix1*) = 
= (— 20x)" XL 
. Cex) er 
2(x1")? Up2 


: "P2 
Poet te 


where 
P< Sa) + 1 
ioe P2 
X1' = 3,487562189 = 
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and X%, P2 must be again determined from program reduction. 





Now 
ey pasa? Mp2, (| __ x1 ie. 
= | 
XL Se See xl 
and 
Uy if 2(X1)* ee 
ee 2 
XD ee = 
_ cn eel 
Eye) 8 
and 
i > 22 
i aon 
where 
. Ove 
Suen 
so that 
U 
Z1 So 1- H H 


In summary, to calculate the uncertainty in Ge at a given point 
the following parameters must be obtained: R2, Ql, V7, X@, 
BY, 28, P2; and Xl'. Note: X1l' must be calculated as its 


value is not retained in the current reduction program. 
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Similarly for the rotor loss coefficient (tp), the prime 
effect is also through the R2 term. We can obtain a satis- 
factory approximation by dropping all contributions to the Ll 


term except X2. This gives 











POMBO eat = (k2)*/,, °x2 

fea [ole aT ee X2 
Therefore, 

eel 

X Ey! 
and 

iy T 

i ae a 

Z So = oa 
where 

. Lo 
> Say 


Using these approximations for calculating the Be and oR 


we obtain for the sample point 


= .593 and = 0996. 


za 0273 
Z1 Z3 
Comparison shows that for Bc this is only 3% below the exact 
uncertainty and for oR only 10% below the exact uncertainty. 
The above approximate formulas have been incorporated into a 


computer program and used to calculate the uncertainty in R2, 


P2, Z1, and 23 for all points in runs 10-15. The program is 
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called "UNCERT" and is listed as part of Appendix A. A 
summary of the approximate formulas and required inputs is 


listed in Table C-IV. 


20k 





A8 
Ag 
BY 
Bl 
B2 
B5 
B6 
B7 
i 
Fg 
Fl 
F2 
BS 
F4 
F5 
Gp 
Lg 


Ll 


M5 


M6 


Ng 


Table C-I 


DEEN ET LONSOER WTERMS 


Partial integral of rotor shroud pressure over 
Partial integral of rotor shroud pressure over 
Equation (8) from Ref. 7 

Equation (16) from Ref. 7 

Beta l (85) - relative flow angle 

Partial integral of rotor shroud pressure over 
Partial integral of rotor shroud pressure over 
Partial integral of rotor shroud pressure over 
iecal axzal force-stator 

Partial integral of rotor shroud pressure over 
Partial integral of rotor shroud pressure over 
Closure Plate force 

Stator exit force 

Stator axial force - force capsule 

Closure plate force ~ force capsule 

Raelo Of “y—-l)y/y = 72857 


Intermediate product in rotor loss coefficient 
Calculation 


Titermediate product in rotor loss coefficient 
calculation 


Rotor torque 
stator torque 
Rotor RPM 


Rotor shroud static pressure 
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ariea 


area 


aked 


area 


area 


abled 


alca 





Nl 
N2 
N3 
N4 
N5 
N6 
N7 
N8 
N9 
Ol 
Pal 


P2 


Po 
Ql 
Q2 
Q3 
Q4 
R2 
R8 
S03 
S4 
S5 
S6 
S7 
S8 


S9 


ReEEOL 
RO GOT 
ROEOr 
RGLOL 
Rotor 
Rete: 
Rotor 
Recor 
ROLOLr 
Rotor 
Ratio 


Ratio 


shroud 
shroud 
shroud 
shroud 
shroud 
shroud 
shroud 
shroud 
shroud 


shroud 


of hood pressure to labrinth pressure 


Pye 
a CS 


TABLE C-I (Cont'd) 


static 
static 
static 
static 
static 
static 
static 
static 
Static 


static 


Average rotor exit 


re 


Racio Of PL Ere 


Hood pressure 


O 


i 


Static pressure at 


Resultant force on 


Static pressure at 


Rotor 


ROCOr 


Rotor 


ROLCOr 


Rotor 


ROCOr 


Rotor 


shroud 
shroud 
shroud 
shroud 
shroud 
shroud 


shroud 


static 
State 
static 
static 
static 
static 


static 


pressure 
pressure 
pressure 
pressure 
pressure 
pressure 
pressure 
pressure 
pressure 


pressure 


pressure 


Stator tap "tip #3" 
Sta vou 
stator tap "hub #3" 
pressure 

pressure 

pressure 

pressure 

pressure 


pressure 


pressure 
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T2 
23 
UL 
U2 
U3 
U4 
Vo 
vl 
v7 
wo 
w4 
xp 
col ' 
Xl 
X2 
x4 
x5 
X7 
X8 
Zp 
Z1 
Z2 


23 


TABLE C-I (Cont‘d) 


qT, 
O 


VO 
Intermediate calculation 
Intermediate calculation 
Dimensionless rotor velocity 
Dimensionless rotor velocity 
Rotor tangential velocity (Vi5q? 
Limiting velocity (V, ) 

Ratio P,/P : 


Xm 


Total mass flow rate 

Ratio W8/g = W9/32.174 

Dimensionless velocity (X54) 

Intermediate calculation for obtaining Xl 
Dimensionless velocity (X4) 
Dimensionless velocity (X,) 

Dimensionless velocity (X55) 
Dimensionless velocity (X_>) 

Blockage factor = l 

Intermediate calculation (y (2(X7) -1)+1) = 2.402 
Intermediate calculation for stator exit pressure 
Stator loss coefficient - - 
Intermediate calculation for X.> 


Rotor loss Coefficient - CR 


Greek letters 


B 


iv 


relative flow angle (degrees) 


ratio of specific heats 
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DABwe C=, (Cont'd) 


C ~ loss coefficient 


8 - ratio of total inlet temperature to standard 
atmospheric temperature (518.7°R) 


Subscripts 

a - axial direction 
ref - referred quantity 
R moo TOtOr 

S =a Stator 

te ae otad lL Conall c1on 
(or t,? 

u - peripheral component 
0 =e Stator entrance 

1 - stator exit 

2 = rotor discharge 


LS 





x 
XL 
X2 
v7 
BY 
Zp 
Bl 
Z2 
X5 
Pg 


P2 


UNCERTAINTIES ASSIGNED 


memperature 


Pressure 


Seatcor Torque 


Rotor Torque 


Force Capsules 


RPM 


0 


Table C=i 


i 


TO RAW DATA PARAMETERS 
eee > 5 
+ 0.04 in Hg 


Gee + 0.54 ain H50 


I+ 


I+ 


RAW DATA FOR RUN 


° 


254.30 


~1442 


Poe 
3571 
1224 
McoG7 05” 
1886 


moO 


OO 29 


17.14645 in Hg 


0 


~ 2427 


Oo br 


Gaeac? Ina br 


= O248elbe 
ee Ome 
Fabken.c=] tr 1 
LOSrPOorNnn 1 

R2e= 130.79 1bt 
i) =1 636.4 °R 
Oly —-— 59.622 in Hg 
M5 = 446.90 in-lbf 
Dome esol. 02 in—ibt 
F4 = -121.87 1lbf 
home) oo0 Lbft 
N = 11977 RPM 
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Parameter 


2 


EZ 


Z1 


Table C-IV 


UNCERTAINTY APPROXIMATION FORMULAS 


Formulas 
UR2 _ 1.9966 
R2 R2 
Ue ed) oR2 
C eee ce 


where C= (BY)? - zp 


Beene (eee HR2, 2 
P2 Bo +SOR(C) R2 
. 5 3elZ Ue, 21/2 
BY +SOQR(C) C 
DME CLG Up? 
D D P2 
2 2 
where D = (X1l')” = (X9)" +1 
ie yoesp ',2 
XI tie er D 
toe x1! eee, 
ly eer a2 
ex) 2 (x1)? xc 
X2 peer) cl 


BINS es (Gah) s 
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Tigh opbye ts 
R2 

BY, 28, 
BY, 29 
BY, 28, 
ancien 
R2a eG 
eb 
Xl*, eg 
Di 
Xf, Xl, 


R2 


Z3 








H _ X2 
a * x3 
2 
where H (X2) Go 
1 - (P2) 
eens nies 2H 
Zl  l-H 4H 
ae : (x2)* Ux 
re pa X2 
ee Se 
Z, eos el 
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ye 


eZee 
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APPENDIX D 


CALCULATION OF CONDENSATION EFFECTS 


IN THE STATOR NOZZLES 


D.1 SUPERCOOLING IF CONDENSATION IS ABSENT 


The present calculations follow the methods given in 


Cahpter 10 of Reference 13. The symbols used (in this 


section only) are as follows: 


td 
Il 


AT = 


Subscripts 
1 = 


> = 


total pressure 


total temperature 


static pressure 
static temperature 
dew point temperature 


Saturation pressure of the water vapor at 
the mixture temperature 


partial pressure of the water vapor 
partial pressure of the dry air 


total pressure of the mixture = P_ + P 


P = wW 
relative humidity = = 
Gg 
Pu 
specific humidity = 0.622 oe 
a 


mach number 


= Top 


compression intake 
turbine plenum 


stator nozzle 


Zo 





The purpose of the calculations which follow is to examine 
the static temperature variation in the flow through the 
stator nozzles and to assess the likelihood that condensation 
will occur. The calculations are carried out for the air- 
water vapor mixture entering the Allis-Chalmers compressor 
at typical Monterey weather conditions. The question is - 
does the normal air flow through the stator nozzle supercool 
sufficiently to produce condensation within the nozzle itself? 
The recorded weather data for Monterey for the days of the 
tests in this report showed a limited variation in relative 
humidity and temperature of » = 60-70% and T = 60-70 °F 
respectively. Therefore, the following conditions were 


assumed: 
6 = 60% = Foe 


The process to be analyzed is shown schematically in 


the following sketch: 





T)=530°R T,=514°R 
P,=14.696 psia P,=29.392 psia 
$=608 
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The calculations are for the worst case, say run 12 or 13, 


where the lowest supply temperature (Ty = 574°R) would give 
O 


rise to the greatest degree of supercooling. 

The calculation is carried out as follows: First the 
pressure and temperature for a perfect gas isentropic expan- 
sion is calculated to an assumed Mach number. Then, assuming 
meme st low remains saturated, but free of condensation the 
partial pressure of water vapor is calculated. From the 
known pressure and partial pressure the dew point temperature 
is obtained from tables. The degree of supersaturation is 
then the difference between the calculated gas temperature 
and the dew point temperature. 

Since no water vapor is added in the compressor, the 
specific humidity of the air does not change from atmospheric 
inlet to stator nozzle conditions. At atmospheric inlet 


conditions, 


Pw 
¢ = 0.60 = ae 
g 
where, at 70 °F, Es = .3632 psia (saturation tables for water). 
Since 
P = 1 atm = 14.696 psia = P_ +P 
W a 
Ps = 14.696 - P 
iL al 
and 
Ee. = 0-60 Ee Seu. oO(. 3632 psa) 


Zonk 





P = .21792 psia. 


ak 
Therefore 
P P 
ee ere G60 a 
oF ‘ P : 14.696 = P_ 
Fit 1 
_ 21792 
= 0.622 (77-496 —31792) 
lbm HO 


POU 26 


lbm air ~ 


This fraction (vy) remains constant up to the point at which 


condensation occurs in the nozzle. Therefore, at station 2, 


y P 
$ ~2 
2 CmacerED 
Jo 
where 
P = 1.4436 at T = 574 °R 
i) 
Since = Yor then 
Pw eae 
pee = on G22 = 
atl 2 
Or 
P P P 
Bee 2 "2 
Po Po ae 
1 9 2 
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Writing, 


Eee 
a 1 7 
Dig a (PS P ) 
eZ. a 2 
li 
wy ey, “wy 
- %. ° *. 
iL 1 
Po aoe Le 
p _ 1 - 
W P 
2 Wa 
e. 
1 


Zope 72/4 696 = 2.21792) 
ee wong 
DAVOS = ee ey 92 


~-43584 psia 


Therefore the relative humidity at stagnation conditions is 


given by 
uD Pa 
ee 009861292392) = 343584) 
2 O56 222 0.622(1.4436) 
32 
= 1, 3019 


Now 9 will increase in the flow through the nozzle from the 
plenum to the nozzle throat to the nozzle exit, since the 
static temperature in the flow progressively drops. 

For the purpose of calculating supercooling in the higher 


mach number region of the expansion, from this point it will 
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be assumed that the flow is fully saturated, that is to say, 
¢ = 1. It will also be assumed that no condensation actually 
occurs, and therefore y = 0.00936. 

The dew points are calculated on the basis of the 


partial pressure, Dis. Therefore, at any station 3 in the 


ae 
mozzie, if = l, 
oe 
So 
oa = 1 
33 
Or 
eh 
3 93 
and 
P P P 
eee. 8 
a a Pe 
ie 2 3 


until condensation occurs. Therefore, 


Eee a 
_ 2 = 2 - 
"Ww 7 Ee i i PS (P, Pw 
Z 2 
or P 
Z 
P_ aa 
Pp = = 
W Pp 
3 Wo 
2 
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Since 


W 
2 _ 43584 _ 
p = 39.392 - .43504 = ‘O15 
a 
2 
.015 P, 
“oe|6oleesos 


when P, 1s obtained from isentropic tables for the expansion, 
the dew point temperature is obtained using Table Bl-b of 
met. 13. 


The results of the calculations for 5 mach numbers (or 


nozzle statics) are listed in Table D-I. 


Table D-I 
CALCULATED SUPERCOOLING (AT) IF NO CONDENSATION OCCURS AT Ty, = 574 °R 
O 
(PL =2 atms, 9=0.6 at intake conditions T=530 °R and p = 1 atms) 
° 
Station Mach Pp oe aT Typ oe 
mee nroat) 41.0 iS. 52774 478 S17 -39 
2 die eZ elie S 7 446 510 -64 
3 3 HeeoO sd 5/4 429 507 -78 
4 (nozzle 1.4 Beez 3G 74 412 S03 -91 
exit) 
5 15 8.006 574 396 499 -103 
[Sample calculation: since P =. 0 lane Pa, acustati1on 2, 
3 
m. peeeor4’ee (12.12) = .1797 psia. From Table Bl=-b of 
3 
Ret. 13, T =O LOR. | 


DE 


Za 


The results of similar calculations for identical condi- 


tions except Ty = 670 °R, are given in Table D-II. 
O 
Table D-II 
CALCULATED SUPERCOOLING OF (AT) IF NO CONDENSATION 
OCCURS AT Ty, = 670 °R 
O 
(Py = 2 atms, >» = 0.6 at intake conditions 
> T = 530 °R and P = 1 atms) 
Station Mach P oe i Typ AT 
1 (throat) 1.0 U5 252 7 670 S56 SL +41 
2 eZ eZee le? 670 SA) 50 +10 
3 es 10.608 670 501 507 - 6 
4 1.4 9236 670 481 503 -22 
5 ies 8.006 670 462 499 -37 


Note that a negative delta T implies supercooling. The 
results for the lower temperature (TY = 574 °R) indicate 
that the flow is supercooled from meine cae to the exit 
plane, with AT = -103 °F at m= 1.5. Although the stator 
nozzle is designed for a Mach number of 1.4, expansion of the 


Flow to the interstage pressure conditions will result ina 


slightly higher Mach number approaching 1.5. 


D.2 EFFECT OF CONDENSATION ON PRESSURE 
The effect of condensation on pressures which would occur 
in the nozzle expansion were calculated using the method given 


on P. 205 of Ref. 15. If the condensation process is treated 
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as one in which the heat of condensation, Beer is added in 
the expansion process and acts to increase the stagnation 


temperature, then the effect on the static pressure is 


given by 
2 K-1 ,,2 e 
Km~ (1 + —— M”) hig: 
(l - M*) pct 


O 


where K = ratio of specific heats and f = mass fraction which 


condenses. Using K = 1.4, 


ee =—elOOOm Btu, lpm; Co = .45 Btu/lbm-°R 
and 


gh = 574 °R , 


Saas Letting M = 1.2, 


dP _ _ (1.4) (1.2)*(1 + 0.2(1.2)*) ,, 1000 
P ery -446 (574) 
= 23.05 £. 


Consequently for a 1% change in the pressure ratio to occur 
dP 
(= = OL 01), 


lbm HO 


A = O4017 23.05 = 0004338 Ibm air 


If this is the fraction of the total gas flow which is required 


to condense, then the fraction of the water present in the 


flow required to condense is given by 
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.0004338 = 
0936 =n, O46336.. 


< {th 


Hence, 4.6% of the water in the flow would have to condense 


at M= 1.2 station to produce a 1% change in the pressure 


ratio. 
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APPENDIX E 


RIG PROCEDURES 


1. Prior to each run, observe and record the calibra- 
mon of: 
a. stator axial force capsule 
b. stator torque capsule 
e.) BOcon torque capsule 
d. closure plate axial force capsule 
e. eee 

2. Ensure aftercooler is drained so that a consistent 
high TT, will be obtained for the run. 

38 Rawsctte and check calibration of the stator axial 
force capsule and BE as before and after each point to 
ensure consistency. 

4. Avoid points in the 17000-18500 RPM range to avoid 
data scatter in this region. 

5. Include smaller increments of KIS, i.e., smaller 
RPM changes between points to more clearly define graphed 
parameters. A spacing of 500 RPM between points would be 
Suitable. 

6. Allow a minimum of 3-5 minutes after recording data 
and setting a new operating point, for conditions to 


stabilize. 
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